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Abstract. Resorting to classical techniques of Riemannian geometry we develop a 
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natural Hamiltonian systems. Then we apply such method to the Anisotropic Kepler 
Problem (AKP) and we prove that it is not analytically integrable. 
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1. Introduction 

Let M be a manifold and let TM and T*M be the tangent and cotangent bundles of 
M, respectively. Moreover let g%j{x) be a Riemannian metric on M. The corresponding 
flow on TM is given by the Lagrangian £ = l/2gij(x)x x 3 . The resulting Hamiltonian 
system on T*M has Hamiltonian H = l/2g lJ (x)piPj. 

Such systems are very important examples of Hamiltonian systems, since many 
problems in mechanics reduce to the study of geodesic flows on Riemann and Finsler 
metrics . For example, with a suitable transformation, the Kepler problem can be 
described as the geodesic flow on a sphere @, §]. 

In this paper we are interested in developing a geometrical method to investigate 
whether systems that can be reduced to geodesic flows on a manifold are integrable. 

We recall that an Hamiltonian system with n degrees of freedom is integrable if 
there exist n first integrals Ii, . . . ,I n with the property that {Jj, Ij} = for all i and j, 
and such that they are functionally independent. In particular the definition above can 
be used when the Hamiltonian describes the geodesic flow on a manifold. 

If M is an analytic manifold and the metric gij is analytic then the flow is said to 
be analytically integrable if the indicated integrals 7i, • • • , I n exist and can be chosen to 
be analytic. 

In U it was proved that the geodesic flow of any analytic metric on an analytic 
compact two dimensional manifold of genus g > 1 is not analytically integrable and 
in || the discussion was generalized to n dimensional manifolds. However it is clear 
that the obstructions to integrability are not only of topological origin. This is because 
different Riemannian metrics can be associated with the same manifold, even if the 
topology (of the manifold) imposes restrictions to the possible metrics. Indeed Donnay 
H and Burns and Gerber [[5] have constructed real analytic metrics on the sphere whose 
geodesic flows are Bernoulli. Knieper and Weiss || constructed real analytic metrics 
with positive curvature and positive topological entropy. Thus, even if the sphere has 
genus g = they found that there are metrics on the sphere such that the geodesic flow 
is not analytically integrable. Hence to refine the criterion of Kozlov (and Taimanov) 
it is necessary to find a geometrical approach that takes into account the peculiar 
metric defined on the manifold. Moreover the topological considerations of Kozlov and 
Taimanov can be applied only to compact manifolds and it would be useful to find a 
technique applicable in the non compact case. The main purpose of this paper is to 
consider such geometrical approach that both takes into account the peculiar metric 
and can be applied to non-compact manifolds. 

The approach we consider lies on the simple idea that if you have an integrable 
system and you add a perturbation that destroys the initial symmetry then, unless a 
new symmetry appears , you obtain a non integrable system. Examples of this situations 
are discussed in || and where suitable generalizations of the Poicare-Melnikov 



method (see for example |TT| , [T2| and references therein) are applied to the perturbed 
systems. In [[| the Anisotropic Manev problem is discussed for weak anisotropics, and 
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resorting to the Melnikov method the occurrence of chaos on the zero energy manifold 



is proven. In || the broken symmetry is the rotational invariance. Similarly in | I0[] , 
where we discussed black holes immersed in uniform electric and magnetic fields, the 
broken symmetry was the rotational invariance. 

To formalize this intuitive idea, in the next section some results of Riemannian 
geometry are recalled, in particular the isometries, the infinitesimal isometries and 
some of their properties are discussed. Then we observe that the equation of Killing 
can be used to find the symmetries that preserve the equation of motion and hence if 
those equation have only the zero solution no first integral can be found other than the 
Hamiltonian. 

In Section 3 we introduce our main example: the Anisotropic Kepler Problem 
(AKP) and we show that with a suitable canonical transformation it can be reduced to 
the geodesic flow on a Riemannian surface with non constant curvature. The Anisotropic 
Kepler Problem describes the motion of one body orbiting around another considered 
fixed at the origin in an anisotropic space (i.e. such that the force acts differently in 
each direction). The AKP was carefully studied in a series of paper by various authors 
(see [13 , |T4}| for an overview of the problem and of the main results), but not much is 
known for weak anisotropies. Indeed even if there was evidence of a chaotic behaviour 
no proof of this fact can be found in the literature (except for the occurrence of chaos 
on the zero energy level that was proved in 0). Moreover the AKP is an example of 
geodesic flow on a non-compact surface and hence it is particularly suitable to illustrate 
some of the peculiarities of the technique under discussion in this work. Hence showing 
the nonintegrability of the AKP for weak anisotropies seems an interesting application. 

The last section is devoted to show that the Equations of Killing do not have any non 
trivial solutions for the AKP and hence that, by the results of Section 3, the Anisotropic 
Kepler Problem has no analytic first integral independent of the Hamiltonian. 

2. Isometries and Infinitesimal Isometries: The Killing's Field 

We will now recall some results and known facts from Riemannian geometry that we 
need in order to prove the main result of the paper; for a more detailed discussion see, 



for example, [15, 16, 17 



Let M be a n-dimensional manifold with a Riemannian metric gij. An isometry of 
M is a transformation of M which leaves the metric invariant. More precisely, given 
the coordinates X\, . . . , x n on M, a transformation Xi = Xi(x~i, . . . ,x n ) is an isometry (or 
motion of the metric) if 

9ij(zi, ...,x n )= gij(xx{x), . . .,x n (x)). (1) 

It is clear that such motions form a group, that is often called group of motion. Moreover 
a vector field is called an infinitesimal isometry (or a Killing vector field) if the local 
1-parameter group of local transformations generated by in a neighbourhood of each 
point of M consists of local isometries. The set of infinitesimal isometries of M, denoted 
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by i(M), forms a Lie algebra. Now a transformation is an infinitesimal motions of M 
into itself, i.e. a transformation that preserves the metric if the Lie derivative vanishes, 
that is, if the equations of Killing 



are satisfied [I7|. In order that there may be a non-null motion we must have the 
additional condition 

QijCe + 0. (3) 

The first fact to verify is that the infinitesimal isometries are all the transformations 
that preserve the variational problem defined by the Lagrangian £ = \j2gijX l x\ so that 
we can reduce the study of the symmetries preserving the equation of motion to the 
study of the isometries. This problem is solved by the following 

Lemma 1 An infinitesimal transformation is an isometry for M with the metric g,ij if 
and only if it preserves the variational problem defined by the Lagrangian £ = \j2gijX l x^ 
on TM. 

PROOF Let us remark that a transformation is an infinitesimal isometry if and only if 
the Lie derivative of the metric is zero, i.e. if the Killing's equations (0) are satisfied. 
On the other hand a transformation preserves the variational problem defined by the 
Lagrangian ( |T3"D when the Lie derivative with respect to the field £ l is zero, 

d£ dP d£ 
that after few manipulations takes the form 

It is trivial to see that (§) implies fl5|). The converse is also true since x l is arbitrary 
and the tensor in the braces is symmetric. This concludes the proof. 

With the preparations above we are well on our way to establishing the following 

Theorem 1 Let M be a analytic n-manifold with Riemannian metric gij and £ = 
\j1gijX % xi be the Lagrangian defined on TM. If the Killing's equations for M have less 
then n — 1 independent solutions the geodesic flow is not analytically integrable. 

To complete the proof of the theorem we can remark that to each integral of motion 
there corresponds a one parameter group of transformation that preserve the equation 



of motion [0 . Hence, since the equation of Killing admit less than n — 1 independent 
solutions, by Lemma [l], there are less than n — 1 transformations that preserve the 
variational problem. The remark above concludes the proof. 
From the previous theorem we easily obtain the following 

Corollary 1 If the only solution of the Killing's equations is £ l = 0, the corresponding 
geodesic flow is not analytically integrable. 
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We now want to consider some more properties of the infinitesimal isometries. A classical 



result |L6], [u| is the following 

Proposition 1 The Lie algebra of infinitesimal isometries i(M) of a connected 
Riemannian manifold M is of dimension at most \n{n + 1), where n = dimM. If 
dim i(M) = ^n(n + 1), then M is a space of constant curvature. 

This Proposition gives the maximum number of infinitesimal isometries and hence also 
the maximum number of first integrals for a geodesic flow on a manifold of dimension 
n. In particular in the two dimensional case we have 

Corollary 2 If M is a two dimensional connected Riemannian manifold of non 
constant curvature then the Lie algebra of infinitesimal isometries i(M) is of dimension 
at most 1. 

One example that illustrates the results above is the Kepler problem on the plane. Such 
problem can be described as the geodesic flow on a sphere S 2 , hence the isometries 
form a Lie algebra isomorphic to SO (3), and there are three first integrals other than 
the Hamiltonian (one component of the angular momentum and two components of the 
Runge-Lenz vector). On the other hand, if we consider, for instance, an homogeneous 
potential, other than the harmonic oscillator and the Kepler one, we can describe the 
motion only as a geodesic flow on a surface that has a 5*0(2) symmetry. 

3. The Anisotropic Problems as Geodesic Flow on a Surface 

The Hamiltonian function for the Anisotropic Kepler Problem on the plane can be 
written as 

#(q> P) = - 7 2 , a 2M/2 ( 6 ) 

where /i > 1 is a constant and q = (91,92) is the position of one body with respect 
to the other, considered fixed at the origin of the coordinate system, and p = (pi,P2) 
is the momentum of the moving particle. The constant measures the strength of the 
anisotropy and for \x = 1 we recover the classical Kepler Problem. 

With some canonical transformations we want to write the problem as the geodesic 
flow on a surface. To accomplish this task we need the following Lemma 

Lemma 2 Let (q(t),p(t)) be a solution of the Hamiltonian system with Hamiltonian 
H(q,p) which lies on the level surface H = 0. Let us change the time variable t to 
r along this trajectory according to the formula dr/dt = 1/G(q(t), p(£)) 7^ 0. Then 
(q( r ); p( r )) — (q(^( t ))? p(^( r ))) is a solution of the Hamiltonian system (with respect 
to the same symplectic structure) with Hamiltonian H = HG. If G = 2(H + a) with 
a = const, one can also put H = (H + a) 2 . 
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The statement above can be found in |I8] and the proof is based on the idea of extended 
phase space (see If we perform the change of time dr/dt = (qf + fxq 2 ) 1 ^ 2 on the 

submanifold H = h, according to Lemma |2], we obtain the new Hamiltonian 

(qf + ^l\E -h) = (qf + ^ (^y^) - a. (7) 

Then again change the time r to s by ds/dr = 2((g 2 + / ug|) 1 / 2 (iJ — h) + a) on the same 
level surface H = h. By the second part of Lemma || we obtain the Hamiltonian system 
with Hamiltonian 

H = fa 2 + /ig 2 2 ) (p2 " 4 2/1)2 • (8) 

Finally, applying Legendre's transformation, regarding p as coordinates and q as the 
canonically conjugate momenta, we can write 

K = f = f(p 2 -2/>) 2 

(9) 

^ = ff = T(p 2 -2/>) 2 

where the prime indicates the derivative with respect to s. Consequently 

- Vi 

= (p2 "" )2 • (10) 

M (p2-2h)2 

Moreover since we have that 

£ = PiQi(p,Pi) +P2Q2(p,p' 2 ) -H(q 1 (p,p' 1 ),q2(p,p' 2 ),Pi,P2) = (H) 
we obtain a natural Lagrangian system with Lagrangian 

L - W^W ( pf + f) (12) 



and if we let pi = x,p 2 = y and x l = (x, y) we obtain 
1 
2 

where 



£ = \gijX l x 3 (13) 




The Lagrangian above defines a Riemann metric of Gaussian curvature given by 

K = -(1- /i)(x 2 -y 2 ) -2h{l + fj) (15) 
that for fi = 1 describes a metric of constant Gaussian curvature (positive for h < and 



negative for h > 0) P, H. 18]. It is interesting to remark that, by Corollary 2, for [i/neql 



we have at most one infinitesimal isometry (and hence also at most one first integral), 
since the curvature is not constant. 
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4. Nonexistence of Solutions of the Killing's Equations And 
Nonintegrability in the Anisotropic Problem 

In the previous section we reduced the Anisotropic Kepler Problem to a geodesic flow 
on a two dimensional surface, and hence we can apply the results obtained in Section 2. 
This section is devoted to prove the nonintegrability of such system. Thus we are well 
on our way to proving the following 

Theorem 2 The Anisotropic Kepler Problem is not analytically integrable, i.e. it does 
not have any analytic first integral independent of the Hamiltonian. 

The only fact that remains to be proved is that the Killing's equations doesn't have 
any non-null solution. The equations of Killing for the metric ( O ) can be written as 

Mi i l§£ = o 

dy (j, dx 

-4 I g^J + 2fi = (16) 

*(x 2 +y 2 -2h) ^ dy U 

Subtracting the third equation from the second we get 

d JL = d JL, (17) 

dx dy 

combining (|17|) with the first equation in system (16) we obtain 



dx 2 dy 2 ' dy 2 /i dx 2 

Now multiplying the second equation in system ([H]) by (x 2 + y 2 — 2h) (assume 
h < 0), differentiating with respect to x and using equation fllTD and the first relation 
in (|l6|) we get the following second order differential equation 

d 2 ? 1 d£ l 

2/i(* 2 + y 2 -2h)^- + 4? = 0. (19) 

Similarly from the third equation in flTBJ ) we get another second order differential 
equation 

2 8 2 f 2 4 df 2 

V +v >-2fc)^- yj- +i e= . (20) 

Lai JU Laj vj Jo 

It is easy to see that if h < 0, y = is an ordinary point since the functions 
—2y/(x 2 + y 2 — 2h) and 2/(/i(x 2 + y 2 — 2h)) are analytical in y at y — 0. Let each 
of these functions be represented by its Taylor series at y = on the real line , then 
every solution of ([19]) on K is analytic at y = and can be represented on K by its 
Taylor series at y = 0. Similarly every solution of ( p0|) can be represented on K by 



its Taylor series at x = since x = is an ordinary point. Therefore we can look for 
solutions of equations fll9[) of the form = X]^=o a n( x )y n an d solutions of (|20f) of the 



form £ 2 = Y2^=o^n{y)x n - Routine calculations show that, from ( ^9|) we must have 

a2 ~~^x 2 -2h) ' " - 2h) (21) 
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and the following recurrence relation 

n - 2 + finjn - 3) 

n+2 ~ ^-2h)(n + 2)(n + lf n [ > 

holds for n > 1. From equation ( pC|) we get 

//6 , (A* - !) & i 



(23) 



(24) 



(y 2 -2h) ' ° 3(y 2 -2h) 
and the recurrence relation 

2fM + n(n — 3) 
"+ 2 ~ _ (3/2_2/i)(n + 2)(n + l) n 
for n > 1. With the preparations above we can find the general solution of (|19[) 

= <*,(*) (l - ^^y 2 + 6fia j^% y 4 + •••)+ ^ 

+a 1 (x) (y - 3i$t=2h)V 3 + 30 ^p-2fc)^ 5 + 
and the general solution of ( |20| ) 

= b (y) (i - + $^x* + ...) + 

We will not need the explicit form of the solutions, but we will retain the recurrence 
relations. Now we need to check that the solutions we found above satisfy the 
"consistency" relations ([181) . Consider a solution of (|19|) written as a power series 
J2^=o a n(x)y 2 , using the first of equations ( pT8[) we find the relations 

d 2 a n 

-—r- = -2fia n+2 (27) 
dx z 

since, by analyticity, we can exchange the sums with the derivatives. It is easy to check 
that choosing n = the previous relation reduces to 



(26) 

3 _l (MzHE ^5 
30(x 2 -2h) 2 X 



d a 2a 



dx 2 ^ 2 [x 2 — 2h) ^ 

that is solved by = (x 2 — 2h). For n = 2 instead we have 
d 2 a,2 „ 1 — u 

From the first of equations ([H]) with do = (x 2 —2h) we get that a 2 = —1 / fi doesn't satisfy 
(|9|) if /i 7^ 1. Thus the consistency conditions impose that ao = and consequently 
&2n = for every integer n. 

Now, to show that the odd terms of the series also vanish, we consider equation 
(f27p with n = 1 and we get 

d 2 a x 2(u - 1) 
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also for n = 3we obtain 

d?ci3 a% , . 

= = 5(x*-2hy (31) 

Differentiating the second of equations (|21|) twice and combining the result with equation 
(|30|) and (pi]) we find the following differential equation 

3(x 2 -2h)^ + 12x^ + 2(2 + M )a 3 = (32) 

Combining (|3l|) and (^) we can eliminate the second order derivatives and we obtain 

the following first order differential equation 

da 3 2(2 + ^) + 3/5 7 /qq\ 

& = Ito ° 3 = ~x as (33) 

where 7 = 2 ( 2+ ^)+ 3 / 5 and the general solution is as = Cx~ n and satisfies equation (|31|) 

if and only if C = 0. Consequently 03 = and hence a^n+i — for ever integer n. 

Therefore £ x = 0. 

Similarly we can apply the same reasoning to prove that £ 2 = 0. Indeed it is enough 
to consider again equations (^-P3|) exchanging x and y, replacing /1 with and a n 
with b n . Therefore we can conclude that there is no non-null motion, i.e. g^^tf = 
consequently, by Corollary 1, the geodesic flow is not analytically integrable and thus 
the initial anisotropic problems does not have any analytic first integrals other than the 
Hamiltonian. This completes the proof of Theorem 2. 
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